  % Load an RF signal sampled at 20 MHz
  load('../RFsignal@20MHz.mat')
  % I/Q demodulation
  [IQ,Fc] = rf2iq(RF,20e6);
  % RF signal and its envelope
  plot(RF), hold on
  plot(abs(IQ),'Linewidth',1.5), hold off
  legend({'RF signal','I/Q amplitude'})
  fprintf("Center frequency: %f\n",Fc);